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1. Introduction 



Finding a consistent and predictive UV completion of gravity is one of the most challenging 
tasks in theoretical high energy physics to date. In this endeavor it soon became apparent that 
understanding the theory's renormalization group (RG) flow may be a crucial ingredient. The em- 
phasis is thereby on fixed points (FPs), which could provide a consistent UV completion of the 
theory within Wilson's formulation of renormalization. In this context, it is natural to distinguish 
the following cases: firstly, there may be a Gaussian fixed point (GFP) linked to the free theory. 
This is the structure underlying the asymptotic freedom of a perturbatively renormalizable the- 
ory. Secondly, the flow may possess non-Gaussian fixed points (NGFPs) where the corresponding 
fundamental action contains interactions. A theory, whose UV completion is provided by such a 
NGFP is termed asymptotically safe. Notably, asymptotic safety may be as predictive as asymptotic 
freedom |jl|]. 

Already at a very early stage, it was observed that the UV-completion of the Einstein-Hilbert 
action is not given by a GFP, i.e., the perturbative quantization procedure does not lead to an 
asymptotically free quantum field theory [^. In order to rescue perturbative renormalizability the 
Einstein-Hilbert action was complemented by introducing fourth order operators This improves 
the UV behavior of the theory and the marginal couplings associated with the four-derivative terms 
are asymptotically free at the one loop level [Q, |5|]. Unfortunately, this improvement comes at the 
price of introducing massive negative norm states ||^ 0], so that it is commonly believed that higher 
derivative gravity is not unitary. 

Along a different line Weinberg proposed that gravity could be asymptotically safe A 
key ingredient in investigating this possibility is the gravitational version of the Wetterich equation 
[^], which allows to investigate non-perturbative properties of the gravitational RG flow. Since 
its advent this tool provided an impressive body of evidence that gravity indeed possesses a suit- 
able NGFP [p^]. Recently, these developments have culminated in the proposal of a systematic 
algorithm for solving the gravitational functional renormalization group equation (FRGE) using 
off-diagonal heat-kernel techniques: the universal RG machine (URGM) [[TH]. 

In the sequel we will rederive the perturbative jS -functions of higher derivative gravity ^ 
from the FRGE [^]. Similar studies have been carried out before [12, 13], where it was observed 
that keeping track of the quadratic and quartic divergences has a drastic effect on the fixed point 
structure of the RG flow. These contributions shift the fixed point for Newton's constant and cos- 
mological constant to non-zero values, rendering the theory asymptotically safe instead of asymp- 
totically free. The main purpose of the present work is the demonstration that the URGM recovers 
these results. Surprisingly, the regularization scheme intrinsic to the URGM unveils certain features 
in the fixed point structure of the theory, that have not been stressed before. 



2. The flow equation for higher-derivative gravity 



In the case of higher-derivative gravity, our ansatz for the average action contains all gravi- 
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tational interaction monomials with four or less powers of momentum: 



IZkAk - ZkR + -^C^ - + 



2Xi 



3A, 



At 



+ ^GF + + 5b . 



(2.1) 



Here = Cf^yajiC^""^ abbreviates the square of the Weyl tensor, E = - IR^yR^"' + ^R^ is 
the integrand of the Euler topological invariant, Z = \/\67iG contains the dimensionful Newton's 
constant G, and all coupling constants are allowed to depend on the RG scale k. 

In order to consistently quantize the theory, we employ the background field method to fix 
the diffeomorphism invariance, splitting the averaged metric g^y = g^y + h^y into a fixed (but 
arbitrary) background metric g^y and fluctuations h^y. The gauge-fixing action ^qf is then taken 
of the form 

Sov = hj d^xy/iF^Y^^Fy, (2.2) 

where F^ = D'^h^y — r\Dyh and the bar denotes covariant derivatives with respect to the background 
metric. Since the gravitational part of the action contains terms with up to four derivatives of the 
fluctuation fields, we also allow for four derivatives in Sgv, employing the minimal gauge [Q] 



b J 



where A = —D^ and the gauge parameters are given by Tj 



1 l+4ffl_ 



l-2£0 



, and V, 



(2.3) 



4 \+(0 ' 3 ' " 'b 

respectively. This choice has the virtue of removing all non-minimal four-derivative terms (as, e.g., 
AD^D^) from the Hessian of the gravitational fluctuations. 

This type of higher-derivative gauge-fixing results in two ghost-terms which take into account 
the Faddeev-Popov determinant [^. The operator F^ leads to a complex pair of ghosts c,c with 
action 



d^xy/gcu A5y 



l-2m 

2(l+£0) 



D^Dy-R\ 



while the contribution of Y^^ is captured by a third (real) ghost field b 



Sb = jj d^^b^.Y^^'by. 



(2.4) 



(2.5) 



The key ingredient for deriving the j8 -functions controlling the scale-dependence of the cou- 
pling constants contained in the ansatz (^) is the FRGE for the gravitational average action [^] 



o>A[<J>,<I>] = ^STr 



(- 



(2.6) 



Here, t = log(^/^o), STr contains a minus sign for Grassmann-valued fields, and = {h^y,c,c,b} 
and <l> denote the collection of fluctuation and background fields, respectively. Moreover, Mk{p^) 
is a (matrix-valued) infrared cutoff which provides a ^-dependent mass term for fluctuations with 
momenta < k^. The interplay between the regulated propagator and the derivative of the regu- 
lator thereby ensures that the trace remains finite for all values of k. In constructing the cutoff we 
follow the URGM and choose in such a way that it provides a mass term to the highest power 



We shall neglect total derivative terms. 
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of the Laplacians appearing in the kinetic terms A i— )• Pk{^) = l^ + Rk{^/k^)^ The profile function 
is taken as Rk{p^) = {k^ - p^)B{k^ - p^). 



Upon substituting the ansatz (|2JJ), the trace in ( |2.6[ ) splits into a gravitational and two ghost 
parts, dtTk = T^^^^ + r= + T'', where 



5h5h 



Tr 



\5b5b 



-1 



-Tr 



5^ 



-1 



(2.7) 



We shall now evaluate these traces employing the off-diagonal heat-kernel methods advocated in 
[ ]r5| , |TT], [l6|]. In this course, we neglect the /:-dependence of all coupling constants inside the traces, 
which corresponds to the one-loop approximation of the flow equation. 

We start with the gravitational trace. Abbreviating pairs of symmetric (external) tensor indices 
with a single label /, j, i.e., hi = h^y, etc. and using the variations [^, the part of the action quadratic 



in the metric fluctuations takes the form^ ^ 8^ [^^f^^ + 5gf] = 
Implementing the prescription of the Type I cutoff detailed above fixes i^f^^ = Kij{Pk{Kf- — A^). 
In order to proceed further, we note that the matrix K is easily inverted. Defining [V^''^)]'^ = 
[/:"1]''d|^^^ and [U]'-j = [K^^fWij, the quadratic fluctuations become 



(2.8) 



Since the "interaction vertices" V and U are of mass-dimension two and four, respectively, the 
inverse of the regulated propagator can be constructed perturbatively in U and V . Neglegting all 
interactions of mass-dimension six or higher, we obtain 



jgrav ^ 



d,Pk 
Pk 



Tr 



d,Pk 



Tr 



+ Tr 



, (2.9) 



where the trace also contains a summation over internal indices. Evaluating the traces utilizing the 



off-diagonal heat-kernel [15, 11] yields 
1 



rpgmv 



(4;r)2 



(2.10) 



which encompasses the well-known result for the (cutoff-independent) four-derivative terms. The 
gravitational trace is found by substituting the properly adjusted and rather lengthy matrices U and 
Vof [|]. 

The inverse (regularized) propagators appearing in the ghost traces and have the form 



(2.11) 



Hn the terminology of this constitutes a cutoff of Type I. 

■^At this stage it is consistent to drop all terms containing derivatives of curvatures, since these do not carry any 



information about the flow of the coupling constants contained in the ansatz ( 2. 1 ) 
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where a is a fixed parameter depending on the couphng co and V^^ is an endomorphism propor- 
tional to the Ricci tensor. The inverse of such an operator can be constructed exphcitly as a power 
series of the curvature and endomorphism tensors by generalizing the techniques [jl^. At zeroth 
order in the curvature, an exphcit computation estabhshes 







-51 - ODaD"—^ . 

Pk ^ ^ Pk{Pk-oA) 



(2.12) 



The curvature corrections to this expression can be found order by order in a systematic bootstrap 
calculation. Including all terms up to &{R^) the result is 

1 



2 J/i 



Pk{Pk-oA) 



+ {[Pk,DaDP] + aD„[D^ A]) {[Pk^DpD"] + aDp[D\A]) 



1 



/f(n-aA)2 



(2.13) 



With this result at hand, it is now straightforward to evaluate a generic ghost trace of the form 
jgh ^ Xj-j^ _|_ V]^^dt^f'. Expanding in V and substituting explicit expressions for the commuta- 
tors, the URGM yields 
1 



(471) 



d'x^ k'{2,-l - ^log(l - a))+k^R (1 + 1,/.- Xiog(i - a)) 



+ 



la- 



log(l 



(2.14) 



where 1//^ = a/(l — a) and V = V^^. Notably, the universal four-derivative piece agrees with 
earlier results on the B4 heat-kernel coefficient of the non-minimal differential operator ( 2.11 ) (with 
Pk{A) = A), providing an independent verification of earlier results [[l7|]. Moreover, one can check 
explicitly, that the a-dependence is such that the limit a — is smooth and the resulting trace 
reduces to the one for the corresponding minimal differential operator. The exphcit expressions 
for and are then obtained from ( 2. 14| ) by taking into account the appropriate pref actors and 
substituting 



1 - 2w) , V^. 



'b — " ' " — 2(1+0))' 

Combining the three traces (p77|), we obtain the final form of the flow equation 



(2.15) 



(4 



133 



-^£ + A(i + 8« + 12«2)/?2 



ZA 



k- 



\2aPl + 12(\-lm)P^ 



k^ 



12(1-20)) 



2^^^(1+200)) , ZA((4+112a))A+<:V5) ^ k^p^ , p 



8o)2 ' 60) 

where (also for later purposes), we abbreviated 

pi=6-96co-A'&ttP- , /72 = 65 + 28a) + 8a)2 

7?4 = 35-218a)-352a)2, p5 = -2-20a), 
P7 = 3 + 26a) - 40a)2 . 



;^2ln (1(1 + «))}/? 

36(ra)-6(T^ln(f(l + «)) 



(2.16) 



p?, = 162-540W, 
P6 = l + 86a) + 40a)2, 



(2.17) 
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3. Perturbative /3 -functions and their fixed points 

The j3 -functions dtgj = ^g. governing the scale-dependence of the coupling constants con- 
tained in ([2.1[) can then be read off by comparing the coefficients of the curvature polynomials 



appearing in (2.16). The j3-functions for the marginal couplings are universal in the sense that they 
do not depend on the regularization scheme. Explicitly, we find 



1 133 3 2 
(iip 10 ^ ' 



1 7(56+1718) 



90 



A, 



(4;r 



(25+1098£O+200ffl2) , 
[T 60 '^^■^> 



Despite our quite different computational approach, this result agrees with earlier computations 
[^, [l^], giving credibility to the off-diagonal heat-kernel methods underlying the URGM. The 
j8 -functions governing the running of Newton's constant and the cosmological constant are most 
conveniently written in terms of the dimensionless quantities G = k^/l6nZk and A = k^-^A: 



AG 



Pg=2G 



6(l-2ffl)^4 (l-2o) 



:P2ln(f(l + w)) 



§SS + mihw(p^ + ~^P2> (1(1 + «)) 



(3.2) 



+ 72i(ra) - AP4) + T92kj (/'S + Ape) • 



These are non-universal in the sense that they depend on the chosen regularization scheme. There- 



fore they are expected to differ from the derivation []12[], which employs the FRGE with a Type III 



cutoff. In particular, the log-terms in ( |3.2| ) are a novel feature in the Type I cutoff computation. 
Their appearance can be traced back to the denominators (/\ — aA) appearing in ( 2.12 ), which are 
absent in the spectrally adjusted case. 

Owed to their key role for studying the renormalization properties of a theory, we close with a 
discussion of the FPs g*, {Pgi\gj=g* = 0), of the j8-functions (3.1) and (3.2). The equation ^xisl) = 
has the sole solution A* = and indicates that Xt vanishes logarithmically at high energies. Thus 
the coupling is asymptotically free. The remaining equations in ( p.l[ ) give rise to the two (well- 
known) fixed point solutions 



FPi,: 



A* =0, 



d* 



-171/56, 



ft) 



1,2 



{-0.00228, -5.47}. 



(3.3) 



Substituting this result into (3.2) we find that only FPi constitutes a FP of the full system 



NGFP 



A* = o, e* 



-171/56, 



ft)" 



-0.00228, A* =0.39, G*=2.39. (3.4) 



This NGFP is UV-attractive in all five couplings. Similarly to the computations []1J], including the 
quadratic and quartic divergences in the computation of the j8 -functions has shifted the couplings 
A*, G* to finite values, rendering the theory asymptotically safe instead of asymptotically free. 
Investigating the fate of FP2, we note that the non-universal j8 -functions are well-defined in the 
region co > —I only. This bound can be traced back to the requirement of positivity of the ghost 



operator (2.11). Since FP2 is not within this bound, it cannot be completed to a FP on the full 
theory space. We take this as a strong indication that this fixed point is unphysical. 
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4. Conclusion and Outlook 

In this paper, we used the exact functional renormalization group equation [Q] in order to 
derive the one-loop RG flow of higher derivative gravity based on a novel resummation method 
inspired by [[l7|]. The main advantage of this technique is that it expresses functional traces of any 
differential operator, in particular higher-derivative and non-minimal ones, in terms of traces build 
from Laplace type operators, that are well known and in general easier to compute. Moreover, it 
allows for a straightforward extension of the one-loop flow to the non-perturbative one. 

Besides the logarithmic singularities seen within dimensional reduction, the functional RG 
scheme employed above also takes into account quadratic and quartic divergences in the regular- 
ization procedure [|l2, 13 1. The latter give rise to a fundamental contribution to the flow of Newtons 



constant and the cosmological constant, whose UV-behavior is then governed by a non-Gaussian 
fixed point instead of the Gaussian fixed point seen within dimensional regularization. Further- 
more, our regularization scheme strongly suggests that, out of the two known perturbative fixed 



points (3.3) for the higher derivative couplings, only one is physically viable. The occurrence of a 
non-Gaussian fixed point also at the perturbative level provides further hints that gravity may be an 
asymptotically safe theory. 



It would be interesting to study the fate of the NGFP ( |3.4[ ) in a more elaborate computation 



which also takes the feedback of the running coupling constants appearing on the right-hand-side 



of the flow equation into account along the lines [18]. We will come back to this point in a future 



publication [19] 
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